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The strong version of the Vafa-Witten theorem is shown may not to hold because the zero con- 
densate from a direct computation of the order parameter is found to be a result on the symmetric 
vacuum. The validity of the Vafa-Witten theorem relies then on its weak version, that the Goldstone 
boson is absent in vector-like gauge theories with vanishing 9 angle. The existence of a charged p 
meson condensate, which violates electromagnetic gauge symmetry, is consistent with this weak 
version of the Vafa-Witten theorem when applied to strong magnetic fields in QCD. 
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The behavior of the vacuum is a very important problem in quantum field theory. The Vafa-Witten theorem [l[ 
' offers a strong constraint on the vacuum of vector-like gauge theories. Recently, It was shown that in a sufficiently 
D | strong magnetic field, the QCD vacuum can undergo a transition to a new phase when charged p mesons condense^]- 
In this phase the vacuum behaves as an anisotropic inhomogencous superconductor that supports superconductivity 
£NJ | along the axis of the magnetic field and generates p vortices. As these initial results are either from an effective theory 
for p mesonsQ or from a phenomenological NJL model Q, an investigation of the underlying fundamental QCD is 
needed. For QCD, the results seem to be controversial: paper [j| claims that in terms of the Vafa-Witten theorem, 

■ charged vector mesons cannot condense out in a magnetic field, whereas a later paper Q states that this is not the 
+f ' case. Indeed, both papers agree that there is no Goldstone boson which originally is taken as a typical signature of 
O r global U(l)i 3 isospin symmetry breaking in QCD for strong magnetic fields. The work Q treats it as indicating an 
D ■ absence of condensation, whereas the work Q explains it as the result of the Higgs mechanism. That is, due to the lock 

i j of the U(l)i 3 with electro-magnetic t/(l) cm gauge symmetry, charged p-mcson condensation will induce spontaneous 
C/(l) em breaking; the photon will then obtain mass by 'eating' a Goldstone boson as its longitudinal component. 
£SJ ' Considering that the absence of a Goldstone boson can have different interpretations, the direct computation of a 
^ . nonzero condensate becomes ever more important. Paper [|| explicitly performs this computation and shows that 
ON ' the condensate is zero, whereas paper @ avoids directly discussing this issue. Nevertheless, nonzero charged p meson 
condensation is preferred in effective theories and phenomenological model computations, and more importantly, 

■ intuition based on the important formula for the p meson mass in an external magnetic field 

m 2 ± (Sext) = m 2 ± - e-Bext (1) 



X 
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, tells us that when the external magnetic field exceeds its critical value m? ± /e, the effective p mass m p ± (i? xt) becomes 
• • ' imaginary. This is a sign of an unstable vacuum similar to the famous situation for the electroweak Higgs potential; 
specifically, the curvature of the effective potential for a charged p meson field near the origin (symmetric vacuum) is 
to be changed from concave to convex if the effective p mass becomes imaginary. As long as the effective potential 
is bounded from below, convex behavior of the effective potential near origin demands that there should then exist 
at least two new non-symmetric vacuums characterized by nonzero condensates. One might guess that, in Q, a 
no-condensate result is due to the fact that there are some unexpected infinities not counted in the computation that 
multiply the obtained zero creating a finite nonzero condensate. Unfortunately, up to now, this kind of infinity has 
not been found. One is then lead to question the reliability of the computation. The purpose of this paper is to show 
that this is really the case, i.e., the isospin violation condensation in QCD with 9 = might not be zero. There exists 
a loop-hole in the direct computation of the condensate either in the original Vafa-Witten's work or in the paper [f| . 

The Vafa-Witten theorem shows that the vector-like global symmetries cannot be spontaneously broken in vector- 
like theories with zero theta angle. The result for continuous global symmetries is supported by two facts: direct 
computation shows a zero condensate for the order parameter, and a combination of various inequalities indicate no 
Goldstone boson. If both of these two facts are true, we call that the strong version of Vafa-Witten theorem holds. 
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As discussed above, the absence of a Goldstone boson cannot be treated as a criterion for the spontaneous isospin 
symmetry breaking in QCD with strong external magnetic fields. In this special case, the validity of the strong version 
of Vafa-Witten theorem relies solely on the direct computation of the condensate in QCD and in the rest of the paper, 
we shall mainly focus on this issue. 

In 2-flavor QCD with an external magnetic field, the charged p meson condensate in terms of quark field ip is 1 
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where r and 7 are the Pauli and gamma matrices defined in the isospin and spinor spaces, respectively, t± = (ti±t 2 )/2 
and 7± = (71 ± ij2)/V%- Subscript '£?' denotes the background external magnetic field, 'tr' signifies the trace for 
isospin and spinor indices, = <9 M — igT a A® — iqA cm is the QCD covariant derivative for the gluon fields A° in 
the presence of an external U(l) cm field A e ™ with g the strong coupling constant, and q = efa + l/3)/2 the electric 
charge for quarks. We take the u and d quarks with the same current mass m. The average in ([2]) is computed for 
the gluon field with the standard QCD path integral measure; ~ is to note that we have ignored the multiplication 
of some irrelevant finite nonzero constant. Because term i] ^ ) 1 _ m only has 1 and T3 components and no t± components 

in the isospin space, then tr[7±r± ^^_ m ] { x i x ) = due to fact that tr/ 3 (r±) = tr/ 3 (r±T3) = 0, where tr/ 3 is the trace 
of the isospin component. This situation is similar to isospin violation discussed in the original Vafa-Witten paper, 
where the condensate is 



(0\iP{x)T 3 i;(x)\0) 
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where = — igT a Ap\ is the pure QCD covariant derivative corresponding to the gluon fields. In ©, the term 
— has no T3 component in isospin space, hence tr[ra — ] (x,x) = due to the fact that tr/ 3 (T3) = 0. To avoid 

trivial zeroes, |f| adds to the Lagrangian an explicit infinitesimal isospin breaking perturbation term lim e _i.o ei^Tip with 
r depending on isospins, spinors and space-time coordinate. Vafa-Witten does a similar thing by separating m u 7^ m<j 
infinitcsimally. Both studies show that the result is stable under infinitesimal isospin-violatcd perturbations. Thus, 
from direct computation giving vanishing condensates, [|[ is at the same level as the original Vafa-Witten discussion. 

To examine whether it is possible to avoid the above no-condenste results for $2$ and ©, notice that in the effective 
theory and the NJL model, the covariant derivative includes the p field as p^r a . This p field dependence of the 
covariant derivative will generate t± and T3 components in the i] ^_ m and - ^ — , thereby prohibiting the appearance 
of a final zero. If one further argues that those effective theories and phenomenological models are not QCD, then we 
can do a similar thing in QCD by applying the technique we developed in Ref.[7| to integrate out quark and gluon 
fields and integrate in bilocal colorless meson fields $(x,x') and II(x,x') exactly, modifying ([2]) to 
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where a and p are isospin and spinor index sets, Tr is the trace over color, isospin, spinor, and space-time indices, and 



pi-p, 



(xi,x' 1; • • • ,x n ,x' n ) is the n-point gluon Green's function. For ([3]), the corresponding result is just to replace 
1± T ± m © by T3 and to ignore the q4- Cra term. We see due to the appearance of isospin field II in i ^ +qj ^ _ m _ n that 
the original argument for a vanishing result now does not hold. 

Suppose (UJ does yield a nonzero result; one could still question whether ^ and © give zero results. Why 
does this no-condensate result, once the path integral is changed to some other form, become nonzero? To answer 
these questions, notice that once spontaneous symmetry breaking occurs, the order parameter should be a multi- 
valued quantity. However, ^ and © themselves only give a single valued zero result for the order parameters 
(0\ip{x)T±'y±ip(x)\0)B and (O| , 0(x)t3'0(x)|O); the expected multi-value result does not show up, while (|4|) does implies 



Here, we perform calculations in Minkowski space, the rotation to Euclidian space being unnecessary for the present discussions. 
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the existence of multi- value result. Here the multivalucdness is associated with possible multiple local minima of 
the effective energy functional. These suggest that J5J and ([3]) are just results pertaining to the symmetric vacuum. 
Without the knowledge of other possible vacuums, it is dangerous to jump to conclusion only with the information on 
the symmetric vacuum. To search for multi-valued results, one needs to find some alternative expression which can 
include in the non-symmetric vacuum effect, like that given by (j4|). The standard way to deal with this situation is to 
introduce an effective action to search for different vacuums and assess which is physical according to its vacuum energy. 
In more detail, suppose <fi(x) is an order parameter such as (0\iP(x)t±j±iIj(x)\0) b or (0\ip(x)T3 , ip(x)\0); introduce an 
external source term J(x)4>(x) into the Lagrangian of the path integral and construct the connected generating 

functional W[J] on J and classical field (j) c {x) = . The condensate is then the value of 4> c (x) when we switch 

off the external source J — 0. This (O|0(x)|O) = (/>c\j_ simply corresponds to the present direct computation of the 
condensate given in $2$ and (j3l). 

In practice, usually we do not know whether W[J] is a single- valued functional on J or not. We usually start compu- 
tations by assuming a single- valued functional, hence 4> c {x) should be single valued. To retain multi- valued condensate 
solutions, instead of directly calculating </> c , one introduces an effective action T[<p c ] = W[J] — J d A x4> c (x)J{x) : and 
expresses the effective action T in terms of the classical field (j> c (x), not the original external source J. It is then 
easy to show g^^j = — J(x), which implies that <fi c (x) is an extremum of T. If J = 0, one only needs to search 
for a real minimum of the effective potential (multiplied by the space-time volume and a minus sign, it is just the 
translational invariant part of the effective action) to obtain the physical condensation. For spontaneous symmetry 
breaking, with the exception of the conventional symmetric solution (j) c (x)\j_ = 0, there usually exist other nonzero 
solutions and these correspond to lower vacuum energies. In this analysis to evaluate physical condensate, the key is 
that the external source should not vanish before obtaining the final physical condensate. For ([2]), the corresponding 
generating functional is 



e w[J±1 = J VA^VipVip exp J d 4 

= J VA c l exp|Trln[0-m+7 ± r±J ± ]+i J d A x [-Jg^G q h} , 
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where G°„ is the gluon field strength and in the second equality we have integrated out the quark fields. The classical 
field p±{x) is defined as 
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Note here the external-source-dcpendent isospin violation term 7±t± J± appearing in the denominator of the integrand 
in addition to the covariant derivative ifi. This makes the result nonzero as long as we do not let the external source 
J± vanish. The effective action is 

T[ P± ] = W[J ± ] - J d 4 x p±{x)J±{x) , (7) 

which satisfies ^j^j = —J±(x). Here we emphasize that once we take J± = at the very beginning in ([5]) and 

as has been done in and we will just get the symmetric solution p±(x)| J± _ = 0. The nonzero solution 

can be obtained only when we come to calculate r[p±], which demands that one must not let the external source 
J± vanish in the process of computation. Instead, the external source must be replaced with its function on the 
classical field J± = J±[p±]- Is this constraint on vanishing external source late so important? What is the effect if 
we replace it with its function on p±l Note that to prove a theorem, we need to investigate all possibilities, while 
to show the theorem is not true, just a counter example is enough. Here we exhibit an explicit example that zero 
condensate from directly calculation of generating functional, once considering its effective action, does lead physical 
nonzero condensations. This just contradicts with the logic of Vafa- Witten theorem that zero condensate from directly 
calculation of generating functional will lead physical zero condensation. To avoid QCD complexities and see clearly 
the core of the problem, we take this counter example as simple as possible, which will describe what happens and 
show strong version of the Vafa- Witten theorem may not to hold. 

We consider a 0-dimcnsion field model with generating functional W[J] given as 

W[J] = 12Ar3 cos(6 + ^7r)[cos(e + \tt) + cos(36)] , (8) 
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where A and r are two real positive numbers of the model, quantity 6 relates these and the external source J through 

36 = arc cos ( — ] . (9) 



irX 



The classical field 6 is defined as 



4> = = 2r i cos(6 + ^tt) , (10) 

where we have used the expression for W[J] given in ([8]) and relation ([9]) to obtain an explicit expression for the 
classical field in terms of 6. Notice that switching off the external source J = corresponds to 



e 



= jr, (ID 
j=o 6 



which allows the classical field to vanish, <p\ Q = 0. This is a single- valued vanishing result, just as in © or ©. To 
generate nonzero (f>\ J _ Q , we first solve (|10l) and express the external source J in terms of the classical field <fi to obtain 

J = 4A(0 3 -3ri0) . (12) 

Noting that as our model is 0-dimcnsional, i.e., there is no space-time integration and the functional differential reverts 
to an ordinary derivative, we introduce the effective action as 

r = W- = 6Ari0 2 - A0 4 , (13) 

where, with the help of (fT2"j) . we can perform some algebraic computations to express T in terms of the classical field 
<j>. The result T given by (|13|) is simply the 0-dimensional Ac/) 4 model with the wrong-sign mass-squared, 12Ar?. The 
vacuum is determined by |^ = and it leads to 

[0 2 -3ri]0 = O. (14) 

<f) = is the solution for the symmetrical vacuum, which is already predicted from (fT0|) and (fTTjl originally by directly 
computing the generating functional ([8|). The other two nonzero solutions <f) = ±\/3r~5 are new ones which cannot 
be obtained from (fT0| and (fTTj) . One can easily check that these two solutions have equal lower vacuum energies 
V = — r|,/ = o in comparison with the zero solution, hence correspond to physical condensates. This toy model thus 
produces nonzero condensates from a single-valued generating functional for the symmetrical vacuum, as long as 
we use the effective action formalism. The reason we call (|8]) and (j9]) the generating functional on the symmetrical 
vacuum is that there exist two other generating functionals on non-symmetric vacuums by changing ^ to 

(15) 

which leads exactly to the same classical field expression (fT0|) and the effective action (|13p . When we switch off the 
external source, J = 0, instead of (fTT|) . these two generating functionals on the non-symmetric vacuums give 

e \j=o = \l > ( 16 ) 

I 2 

which combined with (|10[) just give the nonzero solutions (f> = ±\/3r3. Indeed, the three expressions from (J9j> and 
(|15|) just cover the three different solutions of equation (|12[) . if we want solutions expressing <f> in terms of the external 
source J. <f> as function of J is a three-valued function. In other words, it has three branches: the one given by 
@ corresponds to the symmetric branch, the other two given by (fTSj) to the non-symmetric branches. We refer to 
these results as those on symmetric and non-symmetric vacuums, respectively. Our original generating functional (|8]) 
just chooses the symmetric branch. This toy model indicates that due to the choice of a generating functional on 
the symmetrical vacuum, @ and © produces a no-condensate result does not imply that the physical result must 
vanish. 
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One may criticize the simplicity of this toy model that first it is in O-dimcnsion space and then there is no motion 
for the particle; second there is no path integral and then may be no quantum effect. To overcome these shortcomings, 
we generalize above O-dimcnsional A</> 4 model to D-dimcnsional O(N) A0 4 model by considering following generating 
functional 



,iW D [J] 



[Vfa] cxp 



(x)[-i^ + 6Aft] 



- ^[fa(x)fa(x)} 2 + Ji(x)<j)i(x) 



(17) 



where A and f are two real positive numbers of the model. fa(x) has TV components, Ji(x) is corresponding external 
source, i = 1,2, ... , TV. We discuss the large N limit of the model. Under this limit, there is an overall factor TV 
appears in the exponential of the integrand of the path integral due to summation of different components of <j)(x) 
fields, which will suppress all loop contributions [8| 
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(19) 
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Above equation is difficult to solve, we consider a translational invariant situation that external source Jj and fa_ c are 
all independent of space-time coordinates. Then above equation is reduced to 



J 



4A- 2 
N 4, 



12Af3 



(21) 



where J = \J~J~i~Ti and <fi — U <j>i, c 4>i,c ar e modulus of Ji and <f>i. c in O(N) space respectively. Above equation is the 
same as Eq. (fT2)) if we identify J, </>, A/TV, TVr 2 / 3 with J, (f>, A, r 2 / 3 , respectively. This equation has three solutions if 
we want to express the classical field </> in terms of the external source J. Among them, one solution cannot satisfy 
positivity requirement for J and <f>, the left two solutions are 



2TV3r3 cos(6 + -ir) , 
3 



(22) 
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The translational invariant generating functional of large N limit D-dimcnsional O(N) A</> 4 model becomes 

► oo; Ji(x) = Ji 



W D [J] 



N- 



UNXr S cos(9 + -7r)[cos(6 + -tt) + cos(36)] / d u x 
3 3 



(24) 



Wc see that path integral now generates two different branch generating functionals, one is for symmetric vacuum, 
another is for non-symmetric vacuum. If one chooses the symmetric one as [5| and [lj at the very beginning of the 
computation for condensate, it will lead vanishing result. The difference with [5( and [l| is that now instead of single 
value generating functional, path integral gives multi-value result. The reason present direct computation for this 
D-dimensional O(N) A0 4 model resulting in multi- value result is that path integral involves the degree of freedom of 
order parameter fa, just like While in [j| and [lj], due to choice of fundamental quark as dynamical variable and 
the lack of explicitly extracting out the degree of freedom of order parameter, the generating functional unfortunately 
choose the symmetric branch. 
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In more general situations, if we assume the effective action T[4>] be a single- valued function, then for spontaneous 
symmetry breaking, solving equation = —J{x) to express 4>(x) in terms of J will usually lead to a multi- 

valued function, or a function which has several branches, because spontaneous symmetry breaking demands multi- 
vacuum solutions. Among these branches, one is the symmetric branch corresponding to the generating functional 
on the symmetric vacuum. The other branches lead to generating functionals on non-symmetric vacuums. The no- 
condensate results obtained in [BJ and Vafa-Witten's paper fl[ is due to choosing the generating functional on the 
symmetric vacuum from the very beginning. The proofs in [l|] and Q do not deny the other nonsymmetric vacuum 
(or other branch of the generating functional), while our expression explicitly show the possibility of existence 
of other nonsymmetric vacuum. Then just the discussion on the symmetric vacuum may not represent the correct 
result of the whole system. To avoid that misleading result, the safest way is to calculate the effective action, or 
changing over to the generating functional on the non-symmetric vacuum as done in an effective theory or the NJL 
model or (|6]) , where the path integral is already improved involving degree of freedom relate to order parameter as 
our D-dimensional 0(N) A0 4 model. Taking just an infinitesimal perturbation, as Vafa-Witten and Ref.p had done, 
is not enough to cross-over from the symmetric vacuum to the non-symmetric vacuum. 

With this discussion on direct computation of condensate, we find that the no-condensate result in the Vafa-Witten 
theorem and [j| may not be true due to its prejudiced choice of a generating functional on symmetric vacuum and lack 
the information for other possible nonsymmetric vacuum. Ignoring this possible defect, the Vafa-Witten theorem has 
to reduce to its weaker form, i.e., it only proves that in vector-like gauge theories with # = 0, there is no Goldstone 
boson. The charged p meson condensate in QCD for strong magnetic fields happen to have no Goldstone boson due 
to existence of Higgs mechanism. 

The explanation of generating functional on symmetric vacuum may also help to understand the dispute over the 
parity violation problem in QCD [9[. Once we accept the fact that the generating functional Vafa-Witten used in their 
proof is the one associated with the symmetric vacuum, then the absence of a Goldstone boson obtained from this 
generating functional is also questionable. Considering that up to now, unlike the vanishing condensate case, there is 
no exceptions to the result of no Goldstone boson in various effective field theories and phenomenological models, this 
partial result may still be true. In fact, one can check that original Vafa-Witten's discussion on two-point correlation 
function of isospin current for QCD in a strong magnetic field is still true in presence non-vanishing isospin source 
introduced as in Eq.([5]) and the exponential fall-off behavior for the correlation function remains, thus insuring the 
correctness of the weak version of the Vafa-Witten theorem. 

It should be emphasized that our discussions are not only valid for the Vafa-Witten theorem, but also apply to all 
kinds of direct computations of various physical quantities and all generating functional-based results. That is, once 
the system has spontaneous symmetry breaking or multi- vacuums, one must carefully check the results obtained from 
the direct computation of path integrals to see that these are not performed just on the non-physical vacuum. 
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